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Abstract 

Building on our earlier work and that of Son, we construct string theory duals of non- 
relativistic critical phenomena at finite temperature and density. Concretely, we find black 
hole solutions of type IIB supergravity whose asymptotic geometries realize the Schrodinger 
group as isometrics. We then identify the non-relativistic conformal field theories to which 
they are dual. We analyze the thermodynamics of these black holes, which turn out to 
describe the system at finite temperature and finite density. The strong-coupling result for 
the shear viscosity of the dual non-relativistic field theory saturates the KSS bound. 
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1 Introduction 



The hydrodynamics of cold atoms with interactions at the unitarity hmit (for reviews, see 
[1]) is a subject crying out for an effective strong-couphng description. Experimentally, these 
systems are under extensive and detailed study, and rich data exist. For example, the shear 
viscosity has been extracted [3] from energy-loss during sloshing experiments leading 
to an estimate for the ratio ^ which approaches the bound conjectured by While this 
bound is universal in relativistic systems with classical gravity duals P, [5], the system of 
cold fermionic atoms at unitarity is most certainly not a relativistic one, though it does 
have non-relativistic conformal symmetry (for a systematic discussion, see |i7j). One could 
imagine that the nonrelativistic nature of the system has an important effect on this ratio. 
Indeed, the counterexamples to the ^ bound proposed in [8] arise in nonrelativistic systems. 
Theoretically, however, these systems are hard: perturbative techniques are inadequate, and 
lattice methods are difficult to apply to dynamical questions (though see pj). 

In the relativistic context, an effective strong-coupling description of a CFT can some- 
times be found in terms of a gravitational theory in extra-dimensional spacetimes [9l [10] in 
which the conformal symmetry of the CFT arise as the isometrics of the geometry. Via the 
Hawking phenomenon, the thermal ensemble of such systems is constructed by placing a black 
hole in the extra-dimensional geometry [10] . The rigid structure of black hole spacetimes, 
when combined with the finite-temperature gauge/gravity duality, has led to the observation 
of universal behavior of these strongly coupled gauge theories at finite temperature. For a 
nice review of this work, see |llj . 

In the non-relativistic natural guess for a strong-coupling description is a dual 

geometry whose isometrics reproduce the symmetries of the non-relativistic CFT (NRCFT). 
Such geometries were constructed recently in [111 [13], with the metric taking the fornjl 

= L y-^:^ + -2 ) ■ (1-1) 

Here, x is a vector of d spatial coordinates and z is the dynamical exponent, which takes the 
value z = 2 for the fermions at unitarity. In [121 US] this spacetime was shown to solve the 
equations of motion of Einstein gravity coupled to a gauge field of mass m\ = £ii±^ and 



^For previous appearances of related spacetimes in the pre-gauge/gravity-duality literature, see [141 1151 
[TB] . For studies of the supersymmetrization of the Schrodinger group please see [2T| . See also . 
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a cosmo logical constant A = - — p — and was argued to be dual to an NRCFT at zero 
temperature and zero density. To study via duality an NRCFT at finite temperature and 
density, then, we need to put a black hole inside this geometry. 

In this paper we will construct black holes with the asymptotics of (11.11) . show that 
they arise as solutions of string theory, and identify the specific non-relativistic conformal 
field theories to which they are dual. An analysis of the thermodynamics of these black hole 
spacetimes shows that they describe the dual non-relativistic CFTs at finite temperature and 
finite density, with the previously-studied geometry of [131 112] arising in the zero temperature, 
zero density limit. Along the way we identify a scaling limit which describes the system at 
zero temperature but nonzero density. To produce these solutions, we will use a solution- 
generating technique called the Null Melvin Twist ^331 [351 IMl [371 [Mj, which we will review 
in detail below. This Melvinizing procedure is the sought-for analog of the plane wave limit 
described in the introduction and conclusion of [T3]. 

We should emphasize at this point that the NRCFT describing Lithium atoms tuned 
to a Feschbach resonance probably does not literally have a weakly-coupled gravity dual. 
However, the Lithium system has closely related cousins which do have 't Hooft limits - 
indeed, the NRCFTs dual to our black hole spacetimes are precisely such creatures. Our 
hope is that these ideas will be useful for studying strongly-coupled cold atoms in at least 
the same sense in which the Af = 4 theory has been useful for studying universal properties 
of strongly-coupled relativistic liquids, including those made out of QCD. 

The paper is organized as follows. After reviewing the correspondence proposed in [T2l[T3] 
in section 2, we show in section 3 that it can be embedded in string theory. The solutions are 
constructed using the Null Melvin Twist, a machine which eats supergravity solutions and 
produces new ones. The machine has several dials, which we will tune to various ends. The 
input solution that produces the metric ( 11. ip is the extremal D3-brane in type JIB. When 
(in section 3.3) we feed to the Melvinizing machine the near-extremal D3-brane, we find 
that it produces black brane solutions which asymptote to the spacetimes (11.11) . We then 
provide some rudimentary understanding of the identity of the theory at weak coupling in 
this realization. In section four we analyze the thermodynamics. In section five we compute 
the shear viscosity, and show that the strong-coupling universality found by P [5] extends 
beyond the class of relativistic systems. We conclude with a discussion of interesting open 
questions. The appendices contain the details of the Melvinization process, an argument for 
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frame-independence of the viscosity calculation, and some progress towards a 5d effective 
action. 

2 Schrodinger Spacetimes and Non-relativistic CFTs 

Non-relativistic systems which enjoy conformal invariance in d spatial dimension are governed 
by a symmetry algebra known as the d-dimensional Schrodinger algebra. In addition to the 
usual generators of the Galilei group, i.e. generators of spatial translations, Pj, rotations Mj^, 
Galilean boosts Ki, and time translations, if, the Schrodinger algebra includes a dilatation 
operator, D, and a number operator, A^, whose non-trivial commutators are 

[D,P,]=iPi [D,Ki]=i{l-z)Ki [D,H] = izH [D,N] = i{2 - z)N [Pi.Kj] = -5,,N 

where z, the "dynamical exponent", determines the relative scaling between the time- 
coordinate and the spatial coordinates, [t]=length^. In the special case 2; = 2, the alge- 
bra may be extended by an additional "special conformal" generator, C, whose non-trivial 
commutation relations are 

[D, C] = -2iC [H, C] = -iD. 

In this case z = 2, both D and may be diagonalized, so representations of the Schrodinger 
algebra are in general labeled by two numbers, a dimension A and a "number" C.. For 
fermions at unitarity, this number is precisely the fermion number. 

Motivated by the relativistic AdS/CFT correspondence, it is natural to wonder whether 
there exists gravitational duals for non-relativistic CFTs. By analogy to the relativitstic 
case, we expect such a gravitaional description to realize the symmetry group of the CFT as 
the isometry group of the spacetime. However, since there are now two symmetry generators 
which may be diagonalized and whose eigenvalues label inequivalent representations (in the 
AdS case, there is only one, the dimension), we may expect any spacetime which has the 
Schrodinger algebra linearly realized as its isometry group to be two dimensions higher than 
the CFT, as opposed to one-dimension higher as in the case of AdS. 

Such geometries were explicitly constructed in [13J and [12]. More precisely, these papers 
constructed a d + S-dimensional metric realizing the d-spatial-dimensional Schrodinger group 
as its isometry group, and conjectured the associated gravitational system to be dual to non- 
relativistic CFTs at zero temperature and zero density. The metric they presented appears in 
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fll.ip above. We will refer to these metrics as Sch^_^^, where z labels the dynamical exponent 
and d the number of spatial dimensions (note that d+3 = (d+l)+2); in the special case 
z = 1, Sch]l_^_■^ = AdSd+3- Many of the generators are simple to realize as isometries of this 
geometry. For example, the dilatation D is realized as the simultaneous scaling, 

{t,e,x,r} ^ {A%A^-^^,Af,Ar}, 

while a boost K acts as. 



X X — vt ^ ^ + V ■ X — —t. 
Considerably less obvious, but extremely important to what follows, is the identification, 

N = id^. 

The fact that the number operator in a non-relativistic conformal field theory is gapped (one 
Li atom, two Li atoms, three...) tells us that C, must be periodic. But ^ is a null direction in 
the bulk geometry. As such, we appear to be forced into a discrete light cone quantization 
(DLCQ). This will be made more precise in Section 3. 

At first sight, compactifying ^ may look problematic. For example, this may appear 
to violate boost invariance. However, boost invariance remains unbroken precisely because 
the ^ direction is null; this follows from the commutator [N, Ki] = in the Schrodinger 
algebrE^. Perhaps more troublingly, compactifying ^ would appear to introduce a null conical 
singularity at r ^ oo, which suggests that our metric should not be reliable in the strict IR. 
However, this singularity is unphysical. As we shall see below, the singularity goes away as 
soon as we turn on any finite temperature - the would-be null singularity is lost behind a finite 
horizon which shrouds a garden-variety schwarzschild singularity. Meanwhile, physically, we 
always have some finite T in a realistic cold-atom system, and thus a natural IR regulator. 
Finally, and most sharply, even in the strict T — case, the dynamics will resolve this 
"singularity" in a fashion familiar from the study of null orbifolds of flat space [30 l [3T1 [32l [33] : 
a pulse of stress-energy sent towards large r is steadily blue-shifted until its back-reaction 
is no longer negligible; analysis of the back-reaction then shows that the would-be null- 
singularity turns over into a spacelike singularity shrouded behind a (microscopic) horizon. 
All of which is to say, the strict T — limit of our NRCFT is unstable to thermalization 



^We thank Simon Ross for clarifying this issue. 



4 



upon the introduction of any energy, no matter how smalL The challenge, in both the 
NRCFT and the dual spacetime, is not to crank up a finite temperature, but to drive the 
temperature low. 

Importantly, this metric is not a vacuum solution of the Einstein equations. As a result, 
it is necessary to couple the system to additional background fieldstrengths (a pressureless 
dust and a negative cosmological constant) whose stress tensors cancel the non-zero Einstein 
tensor of the spacetime metricj^l. As we shall see in the next section, this system - a metric 
with Schrodinger isometries supported by background fieldstrengths for massive tensor fields 
- has a natural embedding into string theory. 

Further evidence for the conjectured duality is provided by a comparison of Green func- 
tions for scalar operators as computed in the NRCFT and gravity. The n-point Green's 
function of an operator O in a NRCFT, is determined by the scaling dimension and the 
particle number No [7]. Indeed, as shown by Nishida and Son [7], this is the case in any 
nonrelativistic CFT, since that's what's required to specify a representation of the z = 2 
Schrodinger algebra. The two-point Green's function calculated using the gravity theory [13] 
has the same form as that of [7]. Perhaps not surprisingly, the spectrum of the number op- 
erator in the theories dual to geometries of the form (11. ip is the set of integers, since it arises 
from the tower of Kaluza-Klein momenta in the ^ direction. The name for the ^-momentum 
conservation law in the n-point functions 

h + l2 + ... = l'i + l'2 + ... 

is "Bargmann's superselection rule on the mass" [39]. The stress tensor is an operator 
which commutes with the particle number operator (this is a consequence of the Schrodinger 
algebra). The fiuctuations of the bulk metric dual to the stress tensor therefore have zero 
^-momentum. 

Time reversal is an antiunitary operation, which means that it complex conjugates the 
wavef unctions. In particular, say in a weakly coupled theory with field operator it acts 
on an operator Oi by 

T -.01 = ^ O-i. 

■^The recent papers [191 120] find a solution of the vacuum einstein equations with only a cosmological 
constant - this is just the DLCQ of AdS, with the periodic identification breaking the AdS symmetry 
group to its Schrodinger subgroup. As we will discuss in considerably more detail in the next section, this 
corresponds to a degenerate limit of the backgrounds considered in [13U12j and, more generally, in this paper. 
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This is consistent with the fact that our geometries are preserved only by the combined 
operation 

t^-t, e ^ -e. 



3 Embedding in String Theory 

In this section we will show that the zero-temperature solutions described above have a 
natural embedding into solutions of Type II string theory. For simplicity, we will mostly focus 
on the case z = 2, d = 2, though many other cases may be equally straightforward. These 
solutions may be generated in a number of equivalent ways. One useful technique is the "Null 
Melvin Twist," (which was named in which will be described in detail momentarily. 
We will use this technique to construct solutions which embed the Sch], geometry discussed 
above into string theory, as well as solutions which describe the system at finite temperature 
and chemical potential. A second useful technique, which is in fact completely equivalent for 
the backgrounds we consider, is a simple modification of Discrete Light Cone Quantization; 
this presentation will make the structure of the dual field theory transparent. Let's begin 
with the twist. 



3.1 The Null Melvin Twist 

The Null Melvin Twist is a solution generating technique for IIB supergravity which eats 
known solutions and spits out new solutions with inequivalent asymptotics. Melvinization 
has largely been used to construct gravity duals of non-commutatitve and non-local field 
theories [Ml [351 ESI [371 [3H] ; in the case at hand, the Null twisting produces an extremely 
mild form of non-locality, that of non-relativistic theories with instantaneous interactions. 
Importantly, these backgrounds have the property that all curvature scalars are identical to 
those of the original solution [35j; as a result, the constraints on when the supergravity is 
reliable {e.g. A ^ 1) carry over directly. This will also be clear from our analysis of the 
boundary field theory below. 

Interestingly, some solutions with Sch^ asymptotics have already appeared in the Melviniz- 
ing literature. For example, the T = limit of one of the solutions in [351 [37] corresponds 
to S'c/ig^^, though the form field backgrounds break part of the symmetry groupj. These 



*We thank Mukund Rangamani for pointing out the previous appearance of these spacetimes in string 
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were described as dual to "dipole theories" with a non-trivial star-algebra in the dual field 
theories. Here we will argue that these and some other backgrounds generated by the Null 
Melvin twist are in fact dual to NRCFTs. 

The procedure itself is baroque but elementary. The first step is to choose a IIB back- 
ground with two marked isometrics, which we will call dy and dcf). We then (1) boost along 
dy with boost parameter 7, which generates (in general) a new dydt term in the metric, 
(2) T-dualize a la Buscheio along dy, which generates a new dy A dt term in the NS-NS 
5-field and a non-trivial dilaton profile, gyy j—, Bty and $ ^ $ — ^Ingyy, (3) 

re-diagonalize our isometry generators by shifting dcf) ^ dcf) + ady, which generates a new 
term in the metric of the form ds'^ —*>... + {dcf> + ady)^, then return to our original frame by 
(4) T-dualizing back along dy, which generates dydt terms in the metric and dcf> A dy terms 
in B, and (5) boosting back along dy with boost parameter —7. 

All of this leaves us back in the original frame with a new metric, 5-field and non-trivial 
dilaton, all of which are horrendously complicated functions of the two knobs, 7 and a. The 
final step of the Null twist is to (6) simplify this morass by taking a scaling limit in which 
the boost becomes infinite, 7 cxd, and the twisting goes to zero, a —>■ 0, such that the 
product ^ae^ = (3 remains finite. The result is a new solution of the full IIB equations of 
motion with non-trivial background NS-NS 2-form and deformed metric with asymptotics 
inequivalent to the original solution. 

3.2 Rampaging Melvin Eats Extremal D3-brane, Spits Out Sch\ 

Let's apply this procedure to our canonical example, the extremal D3-brane, a solution of 
IIB supergravity with metric, 

ds^ = ^ {-dr'^ + dx^) + h [dp^ + p^ds%) 
where h = l-\ — x- is the usual D3 harmonic function, and self-dual five-form fiirx 

F^^^ = ^dr Ady A dxi A dx2 A dr + Vt^dO A dcf) A dip A dp A dx, 

with ^5 = I cos 9 cos p sin^ p. We must first choose the two isometry directions, dy and 
dcf), along which to Melvinize. Let's take dy to lie along the worldvolume and dcp along the 



theory. 

^The full Buscher rules, and our conventions for them, are given in an appendix. 
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S^] without loss of generality, we can choose coordinates such that y = x^. A particularly 
convenient (though by no means the only possible) choice for d(j) is given by the Hopf fibration 
^ ^ P2 with metric, 

ds^g5 = (isp2 + {dx + Af 

where x is the local coordinate on the Hopf fibre and A is the 1-form potential^ for the kahler 
form on P^, ie Jp2 = dA. We thus take dip = dx- Note that both dy and d(f) act freely, which 
is important for our solution to remain non-singular. 

Melvin being a very messy eater, we will hide the details of the procedure in the Appendix 
and simply write down the result, which is: 

ds^ = ^{-dr\l+/3^p^) + dy\l-(3^p^) + 2dTdy{(3^p^))+hp\dx + Af 
B = 2pp'^{dx + A) Aidr + dy) 
$ = $0 

Note that nothing has happend to the five-form along the way, since T-duality takes dQ^, 
the top form on the sphere, to dy A dfl^, so that the twist dcp ^ dcj) + [3dy of step (3) acts 
trivially. We thus have in our final solution the same five-form flux as in the beginning, 

F5 = (1 + *)VL^de A d(f) A d'ip A dp A dx 

To locate the Schf hiding inside this solution, a few more steps are helpful. Changing 
coordinates to t = {y + r)/ \pi and i = {y — t)/ -\/2, our background becomes, 

ds^ = ^{p^p^dt^ + 2dtd^) + hp\dx + Af 
B = 2(3p^{dx + A)Adt $ = $0 

Adding back in all the terms we dropped in the first step then gives, 

ds^ = i dt^ + '^dtdi + dx^) + h [dp^ + p^rfsls) 

B = 2(3p^{dx + A) Adt $ = $0 



^We can compute A from the Kahler potential on P^, K — thi'Y^\zi\'^ , ie 



Ai — ir- 



For completeness, and because we had a pointlessly slow search for this data in the literature, we present an 
explicit set of conventions and coordinate systems in an appendix. 
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Finally, taking the near-horizon limit, h —>■ R\/ and switching to the global radial coordi- 
nate r = R\/ p, in terms of which h = the solution becomes. 



R2 



-dt'^ + 2dtd^ + dx^ + dr'^ 



-\- Rj^dsg5 



d2 

B = 2V2A-^{dx + A) Adt $ = $0. 

Upon compactifying on the S^, we precisely recover the Schrodinger geometry with d = 2 
and z = 2, our sought-after Sch^. 

Note, too, the appearance of the parameter /3, which was implicitly set to 1/V2 in the 
earlier results of [121 ES]; by a choice of units. The utility of this parameter is considerable 
in what follows. For now, note that retaining it allows a very revealing limit, i.e. (3 —>■ 0, in 
which the solution above reduces to the extremal D3-brane solution with which we began. 
This suggests that there should be a more intrinsic description of our solution as a garden- 
variety deformation of AdS; we will explore this relation later in this section. 

A note on dimensions. As discussed in Section 2, ensuring the quantization of the spec- 
trum of the Schrodinger number operator N requires compactifying the direction ^, some- 
thing not implemented in the Null Melvin Twist described above; this introduces a new 
dimensionful parameter to the game, the length scale L^. Meanwhile, (p is an angular direc- 
tion along a compact space and so dx is dimensionless, which means a, and thus (3, must 
have dimensions of 1/length. Our solutions would thus appear to have two dimensionful 
parameters, and (3. In the case above, however, the specific values of and (3 can be 
rescaled by rescaling the coordinates as 

t-.(3t, e^r^e, 

leaving only the product invariant. It is this ability to scale awajll P which allowed [T2|[T3] 
to set P = l/\/2. Our extremal solutions are thus parameterized by a single dimensionless 
parameter, Holding fixed while scaling /? to zero gives a particularly trivial background 
which respects the Schrodinger group as its isometrics; we shall return to this example below. 

Relatedly, in the above we have set c to 1; it is easy to reintroduce c by taking dt — cdt. 
Interestingly, the way we are getting a non-relativistic limit is not by taking c — > cxd; rather, 
the asymptotic geometry has an effective Ces ~ which goes to oo as we approach the 



^Scaling away the value of P will not be possible in the finite-temperature solutions described below. 
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boundary. In addition to and (3, then, our non-relativistic theory thus contains a finite 
velocity, c. It is natural to interpret this velocity as the speed of sound, Vg- We will test this 
interpretation in the future. For now, we suppress factors of c, which are easy to restore. 

We have thus constructed a solution of Type IIB string theory of the form Schf x which 
is dual, according to the results of [12], to some theory which respects the Schrodinger 
symmetry algebra, i.e. a non-relativistic conformal field theory with d = 2 and z = 2. We 
note that it is straightforward to repeat the analysis of this section for other IIB back- 
grounds with the requisite isometries, including in particular other-dimensional branes and 
other choices of isometries along which to Melvinize. In the next section, we use the same 
techniques to construct a dual description of such an NRCFT at finite temperate and chem- 
ical potential. 

3.3 Solutions with Finite Temperature and Chemical Potential 

As we saw above, feeding the Melvinizing machine an extremal D3-brane produces a solution 
of IIB string theory with spacetime geometry Schl'"'^ x whose dual field theory is a 
NRCFT at zero temperature. Experience with AdS/CFT suggests that putting the NRCFT 
at finite temperature should correspond to the introduction of a Rindler horizon in the bulk 
spacetime, with modes of the boundary theory thermalized by the Hawking radiation of the 
black hole. So we need to figure out a way to embed a non-extremal black hole inside our 
Schrodinger spacetime. It is natural to guess that feeding Melvin a black D3-brane, which 
shares the asymptotic AdS^ x structure of the extremal D3-brane, should produce a black 
hole spacetime which is asymptotically Schl"""^ x S^. As we shall see, this is correct, with 
one important modification which will become clear after Melvin does his thing. 

Thus motivated, let's Melvinize the black D3-brane solution of IIB string theory. The 
starting solution is, 

ds^ = 1 (-drV + dy"" + dx") + h + p2 [dsl, + (d0 + Af]^ 
where h = is the near-horizon limiio of the usual D3 harmonic function and f = 1 + g = 

4 

1 — ^ is the emblackening factor, together with the usual five-form fieldstrength supporting 



^It is easy to keep the full geometry; we skip directly to the near horizon limit here for simplicity. Keeping 
the 1 leads to a simple modification of the above, a geometry which may be interpreted as a non-extremal 
IIB Fluxbrane; it would be interesting to understand the relation between the NRCFT dual to our geometry 
and the worldvolume theory of the full stringy Fluxbrane. 
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the and providing the 5d cc, 

F(^) = (1 + *)fl5 d6 A d(f) A dij A dfi A dx 



Melvinizing this solution along the Hopf fibre is a straightforward application of the 
recipe described above, so let's jump to the chase and simply write down the final result (for 
completeness, the full computation is presented in an appendix). In String Framqj, using 
the coordinates {t, ^,t} introduced above, the result is 

I / B'^ 1— f \l dr"^ 1 

ds^ = -2f^dt^ + 2didt + -^{dt - dif + ^dx^ + ^ + dsl^ + — {dx + Af . 



r'^K V r"^ 2 ^' J j.2 ' j.2j: ' ' ^ 

(3.1) 

where [dx^ A) and ds^^ are as above and 

4 Q2 2 

/ = 1-^, ir = l + V- 

In contrast to the T = solution and to the AdS black hole, the dilaton now has a non-trivial 
profile, 

$ = --Inir, 
2 

while the Neveu-Schwarz two form takes a slightly different form, 

= ^ ((1 + mt + (1 - m) A {dx + A) 

The five form field strength is again unmodified by the Melvinizing. 

It is a simple but tedious exercise to verify that this is a solution to the full lOD IIB 
supergravity equations of motion. Explicitly, it solves 

p=l,3,5 

where T^^^ is the stress tensor for a minimally-coupled p-form field strength H, 

and T^^^ is the dilaton stress tensor. 

There are many things worth noting about this solution. We focus first on the region near 
the horizon. The component of the gauge field along the null killing vector normal to the 
^The transformation to lOd Einstein Frame multiplies the metric by e~^* = Ki. 
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horizon [i.e. B^-^ = A^-) vanishes at the horizon. This is necessary to have a smooth euchdean 
continuation. The geometry contains a nice Rindler horizon with normal (and tangent) vector 
dr, just as in the pre-Melvin hole - in particular, the would-be null-singularity living near 
r — s> oo arising from compactification of ^ is lost behind the Schwarzschild horizon at r = rj^. 
Near the horizon at r = r^^, it is useful to change coordinates to 



r =: Th — W 



up to an irrelevant constant scale factor the metric looks like 

ds^ = dR^ - R\^dT^ + ^dx^ + ... 

with K = In order for this to produce a smooth euclidean cigar geometry, the euclidean 
time coordinate r must be identified according to mr ~ mr + 27r. This is the same result as 
for the pre-Melvin hole. Before Melvinizing, r was also the asymptotic time coordinate, and 
this gave a Hawking temperature j'V^'^-^'^^''^^ = ^ = -i-. With the Schrodinger asymptotics, 
however, t = ^{t + y) is the natural time coordinate. Different Hamiltonians imply differ- 
ent temperatures. The euclidean continuation of our time coordinate t must be identified 
according to %K^2t ~ iny^t + 2n, and hence the Hawking temperature of our black hole is 

Th = (3.2) 

Note that the horizon is unmodified by the Melvin procedure [36]; only the asymptotics (and 
the relationship to the asymptotic coordinates and horizon coordinates) is changed. 

Next note the relative factor of K between the part of the metric and the {dx + A) 
term: turning on a temperature has squashed the along the Hopf fibre (though, since 
/C varies between 1 and 1 + between boundary and horizon, the squashing is gentle). 
When we compactify to 5 dimensions, then, we should expect a non-trivial profile for the 
scalar field associated with this squashing mode. Note, too, that this squashing breaks the 
isometry group of the sphere from S'[/(4) to a subgroup, and thus breaks the supersymmetry 
of the background accordingl}o. Entertainingly, T-dualizing on this fibre still leaves us with 
a squashed sphere, but the dual dilaton is now constant. This gives a perhaps-simpler HA 
description which may be convenient for various purposes. 



For work on the supersymmetric generalization of the Schrodinger group see |21j 
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This solution also admits a number of illuminating limits. As before, taking (3^0 
effectively un-does the Melvinization, returning us to the ordinary black D3-brane solution 
with an identification. Taking T — > sends this solution to the zero temperature solution 
found above, ie to Sch^ x S^. 

All of which raises the obvious question, what is (3? At T = 0, we saw that (3 could 
be scaled out of any physical question by a rescaling of coordinates, which can be thought 
of as a choice of units in the boundary theory. However, this is not the case at finite 
temperature. (3 thus represents a physical parameter of the finite-temperature black hole 
embedded in asymptotically Schrodinger spacetime. To anticipate what property of the 
spacetime/NRCFT this parameter represent ill], look back at the Melvinization procedure. In 
step (3), a turns on a mixing of the y direction (which will eventually become the ^ direction 
after boosting to the IMF) and angular momentum along the . In finite temperature 
AdS/CFT, angular momentum along the sphere translates into finite chemical potential 
for the conserved R-charge dual to the angular momentum current. Combining the above, 
we should expect the CFT dual to this emblackened Sch'i^"^ to have a finite ^-momentum 
density, aka particle number density, which scales as some power of where the factor 

of r'7} with m > is there to ensure that the density runs to zero as T = with 13 

held fixed since, as we have seen, any finite (3 is unphysical at zero temperature so cannot 
determine any physical quantity like density in the zero temperature limit. 

3.4 In Search of QCP: Finite Density at T = 

The solutions found above have non-zero temperature and density. However, sending T ^ 
appears to send them both to zero. This is bizarre, particularly in a non-relativistic theory 
in which particle-antiparticle annihilation is absent so that the number density should stay 
fixed as we take the temperature to zero. We must be able to find a finite-density zero- 
temperature solution! And indeed it is useful to do so, since refrigeration techniques have 
reached the point that thermal effects on the cold atoms can be neglected for many purposes. 

The answer was already implicit in the last paragraph of the previous subsection: the 

particle number density of our Schrodinger black holes scales as some power of the limit 

we should take, then, involves sending ^ oo (which removes the horizon and sends the 

^^We will verify this through explicit computations of boundary thermodynamic quantities in the next 
section. 
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temperature to zero) while taking /5 — > oo to hold the density fixed. A little experimentation 
suggests that the proper limit is to scale ~ ^ oo, keeping = ^ fixed. 

To define the limit more precisely, we make the following replacements: 

t = V2P^i, e=^^, rH = ^'^'fH, f3 = ^P. (3.3) 

We will show in Section 4 that the particle number density in these units is given by p oc 
Q"^ = -ir, while the rescaled temperature is T = ^(^^-^ . To send T keeping finite p, we 
should take ^ ^ oo holding objects with tildes fixed (we will drop the decorations at the 
end). Note that (13.31) includes the transformation which allowed us, at zero temperature, to 
scale away any finite f3 (notably, in the present limit, (3 oo). 

Ignoring the sphere directions for simplicity, the metric in the scaling limit ^ — oo takes 
a pleasingly simple form, 

d'ho = [-^ + 2dtd^ + ^'r'dej + (3.4) 
while the B field takes the equally entertaining form, 

B = ^(dt + 2fiV^de) A (dx + A) 
and the dilaton remains non-trivial, 

$ = --luK, 
2 

where k = 1 + fi^r^. The five- form, as usual, goes along for the ride. Perhaps unsurpris- 
ingly given its pedigree, but surprisingly given its form (note that as r ^ oo, r^K ~ r^, 
so the ^ direction asymptotes to a finite radius controlled by Q), this background can be 
explicitly shown to solve the full equations of motion of IIB supergravity. We will study the 
thermodynamics of this solution alongside that of the finite-temperature case in the next 
section. 

Note, that in our scaling limit, the horizon has run off to r = oo. Happily, the null- 
singularity observed before near r oo is absent thanks to the r'^dC,'^ term - finite density 
has cut off this singularity. However, the dilaton still grows logarithmically. This means that 
the theory contains a region of strong coupling in the IR part of the geometry, somewhat 
similar to gravity duals of IR-strong gauge theories, such as Dp-branes with p < 3 [40j. 
It would be exciting to interpret this scale-dependence of the string coupling in terms of 
screening in the boundary theory at finite density. 
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3.5 An empty trap: Finite Temperature and Zero Chemical Po- 
tential 

If we can take the temperature to zero holding the chemical potential fixed, it stands to reason 
that we can take the chemical potential to zero holding the temperature fixed. Indeed we 
can! However, something rather stupid happens. To see this, let's run the scaling of the 
last section backward, holding T fixed but scaling 0. It is easy to identify the proper 
scaling here: we must take jS 0. Sending the chemical potential to zero thus sends the 
solution back to the black D3-brane with which we began. Of course, we have compactified 
the light-like ^ direction, so what we really have is a DLCQ of AdS (and hence the breaking 
of the symmetry group from S'0(2,4) to its S'c/ig"^ subgroup). However, as an NRCFT, it 
is rather disappointing - there is nothing in the trap. 

As should by now be clear, the backgrounds we have been studying are intimately con- 
nected to DLCQs. Fleshing out this connection will shed light both on the spacetime solu- 
tions themselves and on the NRCFTs to which they are dual. The remainder of this section 
is thus devoted to an analysis of this connection. 

3.6 The Null Melvin Twist with finite as a modified DLCQ 

The Null Melvin twist has the great advantage of being a concrete tool with which to 
generate new IIB solutions from our tired old examples, and as such has been studied rather 
extensively. However, at intermediate steps the solutions are far from simple, and the physical 
meaning of the procedure is rather opaque: what is the intrinsic relationship between the 
final and initial solutions? 

Happily, there is another way of organizing the argument which is completely equivalent 
and which makes the connection between initial and final solutions manifest, following |34j . 
Let's start by studying the DLCQ of the original solution along the ^ = (y — r)/ \f2 direction 
by requiring all fields $ to be invariant under translation along the light-cone ^ direction, 

where = 9^ is the momentum generator on the light-like direction. This produces the 
solution above at /3 = 0. So: how do we introduce /3? The answer is suggested by step 
(3) of the Melvinization, in which we re-diagonalized our symmetry generators to mix the 
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spatial translation current dy into the angular rotation current d(f) d(f) + ady; this replaces 
momentum along the (boosted) ^/-direction, J„, with the sum of Jy — aJs, where aJ^, = d}_.. 
Boosting back to our original frame and taking the boost large turns this into — ^ae'^J^ — > 
— /5J(^ = Jg. To get the final solution with (3 0, then, we should perform a modified 
DLCQ of the original solution in which we orbifold not by a finite translation by L^J^, but 
by the modified current, i.e. we should shift the light-cone momentum of every field 

by (3 times its charge under the d(j) isometry. We'll refer to this as a DLCQ/3. 

The physical meaning of the generated solution is thus relatively straightforward: our 
Null Melvin Twist, aka the DLCQ/?, is a restriction of the original solution to modes with 
fixed light-cone- and angular- moment J^. Note that the DLCQ/? of zero-temperature A/" = 4 
SYM has a particularly trivial limit in which we hold fixed and scale /5 ^ - this is 
just the usual DLCQ of AdS:^. The fact that this system realizes the Schrodinger group 
as isometries is a direct consequence of the compactification of C,- Unfortunately, this limit 
ensures that the particle density of the groundstate is zero, so, while this does describe 
a NRCFT it describes the system at zero density only. It is thus extremely important to 
preserve the parameter f3 if we want to describe something like "fermions at unitarity" rather 
than "no fermions at unitarity." 

3.7 The DLCQ^ of A/" = 4 SYM 

The virtue of the DLCQ/? prescription for our purposes is that it translates relatively easily 
into the dual field theory. To wit, the NRCFTs dual to our Schrodinger spacetimes arise 
as DLCQ^'s of the boundary theories of the original solution. Importantly, the current by 
which we orbifold is Jp = — (3Jr, where Jr is the U(l) R-current dual to the isometry 
current on the S*^, J^, we used in Melvinization. For example, the NRCFT dual to Schg^^ 
is simply the DLCQ/? of the A/" = 4 SU(N) SYM living on the boundary of the AdSs with 
which we started, with Jr the trace of the cartan of the 5*0(6) R-symmetry (corresponding 
to the Hopf fibration we used in Melvinization). 

This result may also be derived via direct application of the Null Melvin Twist to the 

field theory as follows. Start with the Af = 4 SYM theory and compactify it on a circle 

^^Note that for nonzero f3, Jp is not actually light-like in the bulk. So this is in general a DLCQ only from 
the point of view of the boundary field theory. 
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y c:^ y + 2T(Ly with all fields $ required to satisfy the boundary conditions, 

where is the charge of $ under the specified U{1) subgroup of the R-symmetry group. 
This is equivalent to inserting an R-symmetry- valued wilson line around the compact spatial 
direction [31]. For our special case, we chose the R-symmetry such that all three complex 
scalars in the 6 of SO{Q) carry the same charge; this corresponds to a shift on the Hopf 
circle. Now boost the y direction by 7 ^ oo to make it lightlike while scaling a — such 
that /? = \ae^ remains fixed. At weak 't Hooft coupling, this has the following effect. On the 
potential terms in the Lagrangian it does nothing because the R-symmetry is a symmetry. 
On y-derivative terms it amounts to the replacement 9^$ — * 9^$ — i(3qii^. In terms of ^- 
momentum, the net effect is to shift the moding of each field by a constant piece proportional 
to /?, ie 

27CN 

L^{ii - ipgn) = 27ri £ = __ + 

where N is an integer. This is precisely the DLCQ/3 described above. 

This theory seems remarkably simple, even moreso than the usual un-modified DLCQ of 
J\f=A. To understand why, recall that the usual DLCQ tells us to expand every field in the 
theory in modes along the light-like ^ circle. This leaves a KK tower of massive modes, plus 
a single level of massless modes - the zero modes of - which must be treated with, if not 
respect, at least care. The resulting theory is thus deliciously close to being non-relativistic, 
but the persistence of these zero modes reminds everyone that the theory is really Lorentz 
invariant. To get a truly non-relativistic theory, we would like to lift these zeromodes. But 
since the upshot of the DLCQ/? is to shift the moding of all fields by P times their R-charge, 
that is precisely what the DLCQ/3 does. More precisely, the only zero-modes surviving the 
DLCQ/3 are those of R-scalars, ie of the vector bosons, and these must be dealt with carefully; 
among other things, they generate instantaneous interactions between the remaining non- 
zero modes of the matter fields. The result is a theory with only nonrelativistic excitations, 
with the spectrum gapped by two mass scales, the KK scale and the new scale (3. 

In particular, something dramatic happens when /3 x (^) G Z: $ picks back up a zero 
mode (this is just the fact that the wilson line along the DLCQ circle has phase jSL^). (3 
is thus playing the role of an IR regulator for the DLCQ zero modes generated by a wilson 
line. Another curious feature of this scaling is that in order to excite a single KK mode, we 
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need 0{^^) would-be-zero-modes; this suggests that there is an interesting regime where 
1 ^ (A^) <^ where we can drop the KK modes but we still have a well-regulated theory. 
Understanding the interplay between these two scales in more detail, and especially to see 
it arise in the dual geometry, seems worthwhile; we leave such questions to future work. 



4 Black Hole Thermodynamics in Schrodinger Space- 
times 

4.1 Entropy 

The Bekenstein-Hawking entropy density of the black hole (13.11) is 



Note that the dependence on (3 cancels. To write this in terms of more physical variables, 
we need to relate Ly to L^. What we mean by Ly in the formula (14. ip is the extent of the 
horizon in the y direction when ^ has period L^. To figure out what this is, one need only 
look at the metric near the horizon, and plug in. Near the horizon, the metric takes the 
form: 

ds'^ = dy^ + ... = ^{dt - d^Y + ... 

where ... is terms which vanish when we ask about the invariant distance between two events 
separated only in the ^ direction by an amount dC, = L^. Therefore: 

Using the standard parameter map of AdS/CFT (which commutes with Melvinization) 



and the temperature (13.21) . we have 

8 ^ 

For later comparison, it will be useful to note that in units where IGnG^ = 1, we have 

s = L^TT^T\ (4.2) 
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4.2 A comment on the correspondence for the stress tensor 

We begin with a comment about the mysterious-seeming equation (34) of [12], which built 
on ideas of jisj^ . The last section of [12J contains an assertion about which modes of the 
metric couple to which thermodynamic variables of the boundary theory, which is supported 
by matching to a weakly coupled Lagrangian. This expression can be understood more 
directly as follows. In the standard AdS/CFT examples, fluctuations of the metric which 
have nice equations of motion are the ones with one upper and one lower index. These are 
also the components which couple directly to the boundary stress tensor [291 [25j, ^-e. 

luy = O {h') + / V7 2 (Tbd,):: h-; + {h') 

JdM 

where 7 is the metric on the boundary. Given to determine the perturbation of the 
metric g^u = g^j!) + h^i,, the right thing to do is to symmetrize: 

h^,u = ^ {gf^pK + 9uph''^) ■ 

Adding such a fluctuation to the zero-temperature "schrodinger metric" , and setting 

Ao = hl Ai = hl $ = /i*, B, = hl 

gives Son's equation (34) to linear order in these fluctuations. If /i^ is nonzero, there is a 
nonzero fluctuation of h^^. For general z, and restoring factors of f3, we find to linear order 
in the fluctuations 

ds^ = dsl+(^^ - 2P'^^ dt'+^d^dt+(^^ - dt+^d^+h,,dx'dx^+... 

(4.3) 

So T| is the number density of the field theory. An analogy which is useful for under- 
standing this point is the following. In IIB on AdS^ x S^, considered as a ten-dimensional 
theory with a nine-dimensional boundary, what is the meaning of and T^, components of 
the boundary stress tensor with indices in the sphere directions? The answer is that they give 
R-current densities. This is quite analogous to the statement that r| gives the number den- 
sity, since in our correspondence the particle number density is the density of ^-momentum, 
just as the R-charge density is the density of momentum around the directions. 

Note that the interpretation of %^ and remains mysterious. 



13 



We are grateful to Dominik Nickel and Pavel Kovtun for help in appreciating this equation. 
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4.3 Expectation values of the stress tensor 



Consider any bulk theory where the matter lagrangian doesn't involve derivatives of the 
metric. If the boundary metric is flat, the terms in the on-shell action which are linear in 
the metric fluctuations take the form [29]^ 

hdy = O^h"^) + [ {Q>: - Q6i^) + lMy,ct + 0{h^) 

JdM 

where hdy^ct contains counterterms involving the matter fields. 

is the extrinsic curvature of the boundary, and n'^ is an inward-pointing unit normal vector 
to the boundary. Taking n"^ = —y/g^, the formula to extract the expectation value of the 
field theory stress tensor from the bulk data is 

(Tbdy): = -2V^(ei:-5,^(e + a) + ...) 

where 7 is the metric on the boundary {i.e. the metric on a fixed-r subspace), and a is a 
counterterm coefficient, and .... is the contribution of other counterterms. This quantity 
should have a finite limit as r — > {i.e. as it approaches the boundary). 

In these nonrelativistic systems there is one further complication in the extraction of the 
expectation values of the field theory stress tensor. This is the fact that the description 
of nonrelativistic systems one finds here involves an extra dimension ^, whose momenta 
are associated with the conserved particle number. Since the ordinary stress tensor of the 
nonrelativistic system, which we will denote T, is an operator of particle number zero {i.e. 
it is of the form it is related to the boundary stress tensor which depends on ^ by 

extracting the zeromode. This leads to an extra factor of 

= -2 lim L5V^(ee - 5,^(6 + a) + ...) 

We evaluate the stress tensor expectation values in terms of the five-dimensional descrip- 
tion. It would be a useful check to redo this calculation in ten dimensions. The boundary 
counterterms we include are 

Jhdy 

^■^In what follows we studiously set the bulk coupling /C5 = \^gI 
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where = AaA/sj'^^]^ Because of the asymptotic behavior of the solution, we can replace 
the factors of e"'* with their Taylor expansion about the boundary r = 0; will not matter. 
Note that $ here is a proxy for any 5d scalar quantity which behaves as e~^* = K. 

We find that finite expectation values of the physical components of the stress tensor 
require the addition of some extrinsic terms: 

/ext = / V7nM'^F,^a4e"^*. (4.4) 

J bdy 

Note that to the order at which this term contributes to the stress tensor and free energy, 
we can rewrite 

046"** = 04 + 04^4$. 

This term changes the boundary conditions on the massive gauge field away from purely 
Dirichlet [HI |12]. Using our result in appendix C that the coefficient of the F"^ term in 
the 5d lagrangian is e~t*, we see that the special choice 04 = 1,04 = — | implies Neumann 
boundary conditions on A. Remarkably and mysteriously, it turns out that 04 = 1 is required 
for finiteness of T^, and 04 = — | is required for the first law of thermodynamics to be 
satisfied. 

In a scale-invariant field theory with dynamical exponent z, the energy density and 
pressure in thermal equilibrium are related hy^l 

zS = dV. 

Just like tracelessness of the stress-energy tensor of a relativistic CFT (the special case 
2 = 1), this relation arises as a Ward identity for conservation of the dilatation current. For 
our case with d = z = 2, this implies S = V. We constrain the counterterms to cancel the 
divergences., .-.Oa.d so tKattKeWa.diden..v is satisfiei] . Identifying S = —T/^ and 



We thank Dominik Nickel for pointing out that the A'^ term can contribute. 
^^For the special cases z — 1,2 this is shown in [44j; the formula for general z was derived with Pavel 
Kovtun. 

The conditions on the counterterms we find are: 

ai — —6, ai = — (2a2 — 4a4 — 2) , 04 = 1 
6 

for finiteness, and the Ward identity requires 

= 17 - 18ai - 02 (6 + 10a2) + 12a3. 
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V = %\ we fincQ 

S = V = ^{irTY (1 + . (4.5) 

The numerical factor K depends on counterterm coefficients which are not determined by 
finiteness of the energy, pressure, density, or by the Ward identity for scahng. We will fix 
K below by demanding the ffist law of thermodynamics. N will turn out to be zero, in 
agreement with [l6l H7j . As a small check on our calculation, the action evaluated on the 
black hole solution satisfies 

T (/bulk + IgH + let + /ext) lon-shcU ^ ^-^1-^2, 

as expected for the free energy in the grand canonical ensemble. This equality is true of the 
regulated expressions for any choice of the counterterms. 

Note that the thermodynamic potential densities £,V in a. system with dynamical expo- 
nent z should scale like T'^^^ times some function of the dimensionless ratio ^, in agreement 
with our expressions (14. 5p . In our z = 2 case, the factor of makes up for the dimensions 
of the extra power of temperature. 

As discussed in the previous subsection, the density is determined by {T^)- This gives 

H 

Note that the still- mysterious components of the stress tensor are still divergent. That 
some components of the stress tensor would remain divergent in holographic calculations with 
degenerate boundaries was anticipated in [13]. The fact that the components which are hard 
to renormalize are precisely those whose physical interpretation is unclear is heartening. 

4.4 Comments on chemical potential 

Son [12] showed that the mode Aq of the metric in (14. 3 p is the bulk field associated to the 
boundary number density current. The expansion of the finite-temperature metric (13.11) at 
the boundary gives 



p = 2:^ = -L^{^T)\ (4.6) 

T u Z, 



3r^ 

Comparing with the parametrization of the fluctuations in (14. 3p . we see that Aq = |^-|-(9(r^) 
in our background. This suggests that = ^6^^ determines the chemical potential for the 



^^In these expressions we have divided out a common factor of /C = in all of the one-point functions. 
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number density in this background. To extract more precisely the value of the chemical 
potential indicated by these falloffs of Aq requires a better understanding of the couplings 
of these modes [18]. 

Note added in v2: Following ^61 W7\ . we can use a trick to determine the chemical 
potential which makes precise the comments at the end of section 3.3. The null killing vector 
at the horizon is v oc dr = ^ {dt — d^). If we normalize v so that its component along the 
asymptotic time direction is unity, 

v = dt- d^, 

the temperature of the black hole is given by Th = the surface gravity k is defined as 



This corroborates our earlier result that Tu = Now, the fact that the null killing vector 
at the horizon does not point only in the time direction says that the ensemble to which the 
black hole contributes has a density matrix p = e^''(^^^^); this is the translation operator 
by which the euclidean geometry is identified. In our t, ^ coordinates, this gives 

= -1. (4.7) 

4.5 Comments on the first law of thermodynamics 

The first law of thermodynamics should read 

S + V = Ts + pp. 

Given the entropy density, thermodynamic relations determine pp in a system with these 
symmetries. From the Bekenstein-Hawking formula, we have an entropy density of the form 

s = CiL^T^, 

where Ci is a constant. But the thermodynamic relation s = |^ implies 

V = ^ciL^T^ + po{p) = ^Ts+po{p) 

where the second term is temperature-independent but otherwise thus-far undetermined. 
The scale- invariance Ward identity, zS = dV, then implies 

^ + P= + 1 1 P 
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so 

fip = £ + V-Ts= + -i]Ts + \po 
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For our case d = z = 2, this gives 



PP = 2^iLi:T + -po 



Using the thermodynamic potentials S,V,p extracted from T, the enthalpy (the left 
hand side of the first law £ + V = Ts + pp) is 

S + V = ]^L^{'kT)\1 + ^5^). 

Using (jM]), f|4?71) and (g^), the right hand side is 

Consistency of the first law therefore requires N = 0, and determines the integration constant 
pM = 0. 

4.6 Thermodynamics in physical variables 

By rescaling t ^ t' = at,^ —* C,' = b^, we can change the chemical potential to a value with 
respect to which it is possible to differentiate. In these new coordinates, we have 



ii' - -- T' - f -V - I-i— n' - 2^ — 
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The first law still checks. In order to preserve the dispersion relation 2luj + /c^ = {i.e. 
to preserve the gt^ metric coefficient), we should set a = i. In retrospect, the dispersion 
relation with positive mass should be 2luj = k"^; this can be accomplished by setting instead 
a = this will also make the energy density positive. 

Making the substitution a = —^,b = y/JI, ^ = then, we have 

' ~ A p^ ' ^ ~ 2 /i3 ■ 

A small check on this result is the following. The free energy of a scale-invariant theory 
at finite temperature and chemical potential can be written as 

F = -VTy- 
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The power a is determined by dimensional analysis, and for general z turns out to be 
a = Note that this value implies that z£ = dV, in agreement with the scale-invariance 
Ward identity. Free nonrelativistic gases, both classical and quantum with either statistics, 
in the grand canonical ensemble give a = {d + 2)/2 [44J. For z = 1, p = T'^^^f (^) is the 
familiar scaling {e.g. when — 0). The behavior of a at more general z can be argued as 
follows. With scaling exponent z, temperature (which is an energy), scales with z powers 
of inverse-length. Therefore T^/^ scales with one power of inverse-length. The free energy 
density should scale with d + z powers of inverse-length to make up for the scaling of / dtd'^x. 
This gives a = which agrees with the two familiar cases. 



5 Viscosity 

In this section we will study the shear viscosity rj of the fluid described holographically by 
the metric (13.11) . We will do this using the Kubo formula 

r] = -lim-lmG^{u,k = 0), (5.1) 

1^— >o u 

where Gr is the retarded two point function of the scalar mode of the stress tensor: 

G^(^,^ = 0) = -z j d''xdte'-'e{t){[%y{t,x),%y{0,0)]). 
Here we emphasize that the stress tensor is an operator with particle-number zercLj: 

%,{t,x) = / d^T^,{t,x,0 (5.2) 

It was argued in [T7] that very generally the linearized Einstein equation for (p = hy{u)e~'^^^ 
is the scalar wave equation in the same background. The argument uses only the 50(2) sym- 
metry of rotations in the xy-plane; this symmetry is preserved in our solution. We have also 
explicitly checked this statement using the ten-dimensional IIB supergravity equations of 
motion. 

Note that unlike the familiar case of three spatial dimensions, in our d = 2 example there 
is no third dimension in which to give momentum to hy. However, this momentum must be 
set to zero before taking the u ^ limit in the Kubo formula, and nothing is lost for the 
purposes of studying the viscosity. 

^^We thank Pavel Kovtun for an extremely useful conversation on this point. 
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We will show in the remainder of this section that the familiar relation 

7] 1 

s Air 

also holds in these models. Note that the form of our metric violates the hypotheses of the 
general theorem [6]. It would be interesting to see how much further the assumptions made 
there can be relaxed. 

5.1 Scalar wave equation in the finite-temperature solution 

For convenience, we will discuss this problem in six-dimensional Einstein-frame {i.e. dimen- 
sionally reduce on the constant- volume P^). We show that the answer is frame- independent 
in the appendix. 

The wave equation is 

V 9 

In this metric, 

We will study Fourier modes of the form: 

i.e. we have already set to zero the momentum in the spatial directions and the squashed- 
sphere directions. Note that e, are dimensionless variables, measured in units of the 
temperature (times ■^), i.e. they are the gothic variables of [25l [26] . It will be crucial 
to distinguish C from the variable u conjugate to the asymptotic time coordinate t. 

The wave equation becomes 

= u'd^ [^djK^ - (-ye^ + ^'(e - %r + fK. 

The indicial equation near the horizon arises from setting fx = (1 — u)" and demanding 
that the most singular terms at m = 1 cancel. This gives 

e2 

= + -. 

4 
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The solution obeying incoming-wave boundary conditions at the horizon takes the form 

where Fk is analytic at m = 1. Next, to study the hydrodynamic limit, we can expand Fx 
in a small-frequency expansion: 

fM = (1 - (1 + cFi(w) + q,F2(w) + ...); 

here the ellipses denote terms of order C^,Cqy,q^o Plugging back into the wave equation, 
we find, just as in the AdS black hole p8l[24l[25] . 

F,iu) =zln^, F,{u) = Q. 
Using (yf"" = and = , this produces a flux factor 

urjj \41 — -u 4:1 + u J ^ ^ 

where 



167r2 

is the normalization of the bulk action, written here in terms of field theory variables. It 
will cancel mr]/ s. We need the relationship between the momenta associated to the horizon 
coordinates and asymptotic coordinates: 

^"^ " 7f ''^ ^ 

Note that we have restored factors of ^— in the definition of the t, £ momenta relative to 
the gothic momenta. At the boundary -u = e, the flux factor is therefore 



We dropped contact terms in this expression. The real-time AdS/CFT prescription of [2 
says that the retarded Green's function is obtained from the flux factor by 



G^{uj,k = 0) = -2T\u= 



Actually, the correct expansion treats t as the same order as q^; this will not affect the viscosity calcu- 
lation 
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At this point, we pause to consider whose Green's functions we are studying. The 
momentum-space correlator in the Kubo formula flS.ip has had a factor of the volume of 
spacetime divided out by translation invariance: 

the factor VT = L1L2T is 6'^^^{0) in momentum space. As emphasized in equation fl5.2p . 
the field theory stress tensor is the zeromode in the ^ direction of the operator to which /i^^ 
couples. Therefore, when we relate the two-point function of T{t,x,C,) to the momentum- 
space Green's function G^, we should not divide out by the associated factor of L^: 

{[%y{u;,k),%yi-u;,-k)]) = j diidi2 j ti"ftite"^-^^([T,,(t,f,6),7;,(0,0,e2)]). 

Putting this together, the Kubo formula for the viscosity then gives 

77 = - lim -ImG^iuj, = 0) = = 

Note that in d spatial dimensions rj indeed has mass dimension d. This is identical to the 
familiar A/" = 4 answer except for a) the interpretation as the viscosity of a theory in two 
spatial dimensions, and b ) the factors of 1/ \/2 which come from the relation between the 
asymptotic time coordinate and the coordinate which becomes null at the horizon. 

Taking the ratio ^ reproduces the KSS value 

7] 1 

S 4:71 

6 Discussion 

Our black hole lives in a space with very different asymptotics from AdS. There structure 
of the horizon, however, is the same as that of the AdS black hole; this is guaranteed 
by the manner in which it was constructed [36j. The calculation of the viscosity is not 
obviously determined only by the geometry near the horizon. However, the factors conspire 
mysteriously to preserve the viscosity ratio. Our result, then, is some further indication that 
the membrane paradigm should be taken seriously. 

In this paper we have focussed on an example with dynamical exponent z = 2 in d = 2 
dimensions, which is related to the J\f = 4 theory by a twisted version of discrete light cone 
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quantization. Work on constructing string theory realizations for critical phenomena with 
other values of z, d is in progress. 

It would be interesting to find the black hole solution which asymptotes to the NR metric 
with spherical spatial section, i.e. the analog of the black hole in global coordinates in AdS. 
The melvinization can't work quite the same if the starting point is AdS in global coordinates, 
because the analog of y is then an angular variable. 

Having identified a zero-temperature background with nonzero density, and its likely 
weak-coupling description, we can calculate the Bertsch parameter (see e.g. for this 
theory. The Bertsch parameter is the cold-atoms analog of the famous |-ratio of strong and 
weak coupling free energies in the J\f = 4 theory. 

The boundary field theory we are studying clearly contains bosonic excitations, which 
carry charge under the number-density operator. There should be a chemical potential to 
temperature ratio above which they simply Bose condense. In this regard, it would be 
interesting to Melvinize the Sakai-Sugimoto model [IS]; it has a better chance of describing 
a system conaining only fermionic atoms. 

A nice check on our result for the viscosity and our understanding of the thermodynamics 
of the solution will be the location of the diffusion pole in the shear channel of the stress- 
tensor correlators [49] . 

Note added: When this work was substantially complete, we learned that two other groups 
[1^ W7\ had found results which overlap with ours. 
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A Details of Melvinization 

In this appendix, we review the Null Melvin Twist, as formalized in a seven-step dance in 
[36]. 

A.l Buscher Rules and Conventions 



/ _ 1 n' — ^"^y n' — n 9ay9yb + BgyByh 

9yy-— ^"S/ ~ 9ab-9ab - U-iJ 

yyy iJyy 



9yy 



A. 2 The Hopf Vector on 



In constructing our solutions we were forced to pick an isometry direction along . A 
particularly convenient choice involved realizing as a Hopf fibration over P^, which we 
now review to make your life easier than ours was (if you don't already know this stuff). 

The round metrics on P" and S"^""*"^ may be elegantly expressed in terms of the left- 
invariant one- forms of SU(n). For SU(3), these can be written in coordinates as, 

cTi = -{([9cos{ip) +d0sin(6') sin{ip)) 

(T2 = 2('^^^^^('^) ~ d0cos(-?/') sin(6')) 

£73 = ^(d7/; + d0cos(0)) 
In terms of these 1-forms, the metrics on P^ and may be written, 

(isp2 = djJ^ + sin^(/i) [al + 0-2 + cos^(yu)cr|) 

ds% = dsp2 + (dx + sin^(/i)(T3) ^ 

where x is the local coordinate on the Hopf fibre and A = sin^(/i)cr3 = ^^^^-^ {dtp + dcp cos{9)) 
is the 1-form potential for the kahler form on P^ {dx + ^ is the vertical one-form along the 
Hopf fibration). This explicit coordinate presentation is necessary to verify that our various 
solutions in fact solve the full lOd IIB supergravity equations of motion, and to study the 
linearized equations of motion for the fluctuations. 
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A. 3 Constructing the finite temperature solution 



We now walk through the melvinization of the black D3-brane in all its majesty. 
Step 1: We start with the black D3-brane solution, 

ds' = i {-dr'f + dy' + dx') + h(^^ + p' [dsl. + {d<p + Af] ) 

where — 1 + is the usual D3 harmonic function and / = = l — ^is the 

emblackening factor. In what follows, nothing untoward will happen to the dx^, dp^ or ds'^2 
factors, so we'll drop those terms and reintroduce them after the dust settles. The truncated 
metric is thus, 

ds' = i {-dr^f + dy^) + hp\d(f) + 

f t 

Step 2: Boost by 7, ie r — > cr — sy with c = cosh(7) and — s'^ — 1: 

ds^ = i {-dT'^{l + gc^) + dy'^{l - gs^) + 2dTdy{gcs)) + hp^{dx + Af 

f V 

Step 3: T-dualize along the dy isometry using the Buscher rules listed above: 

„2 j_2 / 



ds' 



-dr' 



B ^2dyA dr 



h{l-gs^) 
-gcs 



,2/ 7 , , /i\2 , ,2 1 



+ h[ p\dx + A)+dy' 



$ = $0 - - In 



1-gs^ 

h 



:i.2i 



Step 4: Shift the local 1-form dx to dx + otdy to give 



ds" = -dr' 



f 



h{l-gs^) 



h(p'{dx + Af + dy' 



1 + p^a^{l - gs^ 
1 - gs'^ 



+ 2dy{dx + A){ap'') 



Note that a has dimensions of 



length ■ 



Step 5: T-duahzing back along dy gives 

dr^ 



ds' = 



B 



h{l - gs^ 
dy' 



f 



g^c^s'^ 



+ 



1 + p'^a^il — gs^ 
I- gs"^ 



+ 



2dydT 
h 



gcs 



h 
ap" 



1 + — gs^) 



+ hp\dx + Ay 



1 + p^a^{l - gs^) 
1 

1 + p^a^(l — gs'^ 



1 + p^q;^(1 — gs^) 



{dx + A)A [gcs dT+{l- gs')dy] 



$0 In 

° 2 



1 - gs" 
h 
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Step 6&7: We now boost back by —7 and take a double scaling limit a —> with ac = (3 
held fixed. Since many terms do not survive this, it is easiest to do both steps at once and 
report only the result, adding back in all the terms we dropped in the first step, 

1 



B 



hK 



[-dT\l + /3V)/ + dy\l - f3'p'f) + 2dTdy{(3'p'f)] 



+ —dx'^ + h 
h 



d^ 
f 



+ p^dsl 



Pi 
K 



{dx + Af 



2Pl 
K 

$n - 



{dx + A)A{fdT + dy) 



InK 



Note that jS has dimensions of 



length ■ 



Step 8: Finally, we take the near-horizon limit, h — ^ R\/p^. To compare with the 
solutions of [121 [131, it is convenient to switch variables to the radial radial coordinate 



r 

Ra 



Ra 
P 



in terms of which the boundary is at r = and the horizon at r// = R\/Rh- In terms of r 
and the parameter A = (3R\ we have 



A 2 ^2 ^4 



2^2 



K 



H 



with the metric taking the form. 



ds' 



A 



A' 



,A' 



dr'il + —)f + dy\l - ^/) + 2dTdy{-^f) 



B 



+Kdx^ + K— + {Kdsl^ + {dx + Af) 
2A^{dx + A) A if dr + dy) 



-InK. 
2 



Between the boundary and the horizon, K varies smoothly between 1 and 1 + Impor- 
tantly, the surface r = rn, where / ^ and Bt 0, remains a non-singular null horizon. 
Near the horizon, dr is a timelike killing vector which is perpendicular to the null geodesies 
which span the horizon. We thus have a non-rotating black hole with = 0. This might 
seem somewhat miraculous, since the geometry is not static but, like Kerr, only stationary, 
and so we might reasonably expect a Killing horizon outside the black hole. In fact, this 
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construction, which preserved the near-horizon geometry at each step, had built into it that 
the horizon would be irrotational (and, in particular, have no additional killing horizon). We 
could introduce rotation by starting with a bifurcate killing horizon surrounding an ergo- 

sphcre - i.e. by starting with a rotating black D3 - but, since we will exploit the unbroken 
rotational symmetry of our solution to compute the viscosity, we'll leave this generalization 
to later consideration. 

The upshot is that we have a two-parameter family of finite-temperature solutions labeled 
by the th and A defined in units of Ra- This family has two simple and familiar limits, 
A — > and r^j — > oo. Taking A ^ 0, which sends X — > 1, is easily seen to return us to the 
non-extremal black D3-brane solution with which we began. 

Taking r^^ — oo, by contrast, takes us to the globally non-singular Schrodingcr geometry. 
To see this directly, it is useful to work in light-cone coordinates t = {y — T)j\f2 and 
^— (?/ + ''')/ V^, in terms of which the solution becomes. 



f?2 



■^fdf + 2dtdi - ^{dt - dif + Kdx^ + K— + {Kdsl, + {dx + Af) 
r 1 J 



B = ^A%{dx + A)^{[l + f]dt+[l-f]di) $ = $o-^lnii'. 

In the limit r^j — > oo, which takes / ^ 1 and X — > 1, the metric reduces to, 

2A2 



F?2 



-dr + 2dtdi + dx^ + dr^ 



-\- Rj^ds 



d2 

B = 2V2A^{dx + A) Adt $ = $0, 



which, upon compactifying on the S^, is the Schrodinger geometry with z = 2,d = 2. 
Studying the finite-rij solution near r <^ rn gives the same result. We have thus embedded 
a black hole in an asymptotically Schrodinger spacetime. 

One final set of coordinates will be useful in the computations below. In terms of the 
dimensionless quantities u — r'^/r'jf = Rj^/p'^, 5 = A/th = /3Rh and p = Ra/tr = Rh/ Ra-, 
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the solution takes the form, 



/ di' + 2dtd^ -^(dt- dO'' + Kdx' 

uK \_ u 2 

+ # {Kdsl. + {dx + Af) 
B = y/2S^{dx + A)A{{l + f)dt+{l-m) 



where 



/ = 1 - -u^ K^l + S\, 
These variables simplify many of the computations. 

B Frame (in)dependence of the viscosity calculation 

After compactifying to D dimensions, the string frame metric is related to the I?-dimensional 
Einstein-frame metric by the Weyl rescaling 



In our solution, the dilaton is 



so we have 



_ J_ 

K 



E,D^j^^ (sir) 



In the special case D — 10, this says 5^^,;^° = Klgj^^'^\ 

Now consider the wave equation in a conformal frame reached by an arbitrary power of 
K, where the metric is: 

We have 

jy-5a-l 

det g'^ = det g^^''^^ , ^ = ^^vo1io-d 



where vo1io-d is the constant volume of the compact dimensions, which will scale out of the 
wave equation. 
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The wave equation for a scalar in this background is 

1 / l{aD/2-\ 

^ fK<'^-^y^f 



J^aD /2-1 " 



du du(f) + ... 



U 



The einstein-frame condition above says that in /^-dimensional Einstein frame, K'^ = e ^-^ 
which says 



So we see that in einstein frame, in whatever number of dimensions we want to live in, say 
10 or 6, the factor K does not appear in the wave equation. 

This in turn implies that the viscosity is independent of 5 = 



C Comments on reduction to five dimensions 

Let r = —\\n.K] this is the profile for both the lOd dilaton and the KK scalar associated 
to the Hopf direction. The following two equations are true: 

= lea^ (e(3-i)r^^A^-a,r) + ^ {e^'^-^^^F^ + 2z{z + d)e^'"'-'^^A') 



where 



p2r = J_ 
- K 



ifdr + dy) 



and 



These are the respective equations of motion for A^, and $ for a five-dimensional action of 
the form 
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with m\ = ^^^"^ as usual, and ai + 61 = 03 + 63 = 3z/ — 1, 02 = — 3 and 02 = 03. ^ = | 
is 5d Einstein frame. Here a is the other scalar arising from the KK reduction. The . . . 
indicate terms that do not depend on $, A. We have not yet been able to determine the rest 
of the action. 
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